Superfluid transport of information in turning flocks of starlings 
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Collective decision-making in biological systems requires all individuals in the group to go through 
a behavioural change of state. During this transition, the efficiency of information transport is a key 
factor to prevent cohesion loss and preserve robustness. The precise mechanism by which natural 
groups achieve such efficiency, though, is currently not fully understood. Here, we present an exper- 
imental study of starling flocks performing collective turns in the field. We find that the information 
to change direction propagates across the flock linearly in time with negligible attenuation, hence 
keeping group decoherence to a minimum. This result contrasts with current theories of collective 
motion, which predict a slower and dissipative transport of directional information. We propose a 
novel theory whose cornerstone is the existence of a conserved spin current generated by the gauge 
symmetry of the system. The theory turns out to be mathematically identical to that of superfluid 
transport in liquid helium and it explains the dissipationless propagating mode observed in turning 
flocks. Superfluidity also provides a quantitative expression for the speed of propagation of the 
information, according to which transport must be swifter the stronger the group's orientational 
order. This prediction is verified by the data. We argue that the link between strong order and 
efficient decision-making required by superfluidity may be the adaptive drive for the high degree of 
behavioural polarization observed in many living groups. The mathematical equivalence between 
superfluid liquids and turning flocks is a compelling demonstration of the far-reaching consequences 
of symmetry and conservation laws across different natural systems. 



Consider a flock of starlings under direct threat from 
a peregrine falcon. To dodge the attack, every split sec- 
ond the group collectively takes the decision to change 
direction of motion pQ. Each such change, however, puts 
the flock in a vulnerable condition, which the predator 
is precisely there to exploit. The slightest uncertainty 
may decrease cohesion, or even split the group and push 
some birds astray, leaving them easy prey of the falcon, 
thus decreasing the fitness of the group. This is a gen- 
eral issue for collective decision-making in social species, 
for which forming cohesive groups is a matter of fitness, 
either for anti-predatory reasons, or other environmen- 
tal concerns [2]-[5] . Irrespective of what is the consensus- 
forming mechanism leading to the decision, its actual ex- 
ecution cannot be instantaneous, as a certain time lag is 
needed to propagate the decision throughout the group. 
During this time there is a transient mixture of indi- 
viduals who have already changed state and individuals 
who have not yet done so. For this reason a behavioural 
change of state is intrinsically in conflict with cohesion. 
Consensus must be tight and the decision must spread 
across the group swiftly enough to guarantee robustness 
of the group [7J [8] . 

A collective change of state may be the result of a 
perturbation hitting most individuals in the group. For 
example, a shot heard by an entire flock of birds sitting 
on a tree makes them all take off at the same time. When 
this happens there is hardly any transport of informa- 
tion. More interesting is the case when the collective de- 



cision has a localized spatial origin, starting from a few 
individuals close to each other. This may be due either 
to an external stimulus (a predator, for example) or to 
some spontaneous behavioural fluctuation. In this case, 
the information to change state must spread across the 
group and reach all individuals. It becomes therefore 
essential to understand the mechanism of transport of 
the information. Does the signal get attenuated in space 
and time? What are the laws of propagation and how 
do they depend on the parameters of the group? These 
questions have a major impact on the functional effec- 
tiveness of collective decision-making. However, little is 
known about these problems, both from the empirical 
and the theoretical point of view [9HTT] . 

Here, we perform an experimental study of collective 
turns in natural flocks of starlings (Sturnus vulgaris). 
Studying animal groups performing global turns is natu- 
ral for several reasons. First, turns are a paradigmatic 
example of collective change of state, often taking place 
under severe environmental constraints. Second, turns 
are relevant from an adaptive perspective, as the mech- 
anisms regulating global change of direction have an 
important anti-predatory value in many social species. 
Third, collective turns can be defined sharply from a be- 
havioural point of view and they are relatively easy to 
study empirically in natural conditions. 

We find that collective turns in starling flocks start 
from a few individuals and then propagate to the rest 
of the group like undamped sound waves, a phenomenon 
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FIG. 1. Birds trajectories and turning delays, a, Reconstructed 3d trajectories of three birds belonging to a flock 
performing a collective turn. Sampling at 170Hz we capture fine details of the birds movement, such as the zig-zag due to wing 
flapping (10Hz in starlings - inset), b, c, Trajectories of all N — 176 birds of the same flock as in panel a. Each trajectory 
lies approximately on a plane, justifying a simplified planar description of the velocity, d, The radial acceleration of a turning 
bird displays a maximum as a function of time. In principle, given two birds i and j, one could simply define the turning delay 
nj as the time shift between the peaks of their accelerations. In practice, due to experimental noise, using just one time point 
(the peak) gives an unstable estimate. To calculate Tij in a robust way we must use the entire trajectories. This can be done 
by asking what is the delay Tij by which we have to time-shift the radial accelerations dj(i) to maximally overlap it with di(t). 
This optimal shift corresponds to the time where the correlation function Gij(r) = J dt a<(t) • a,-(i — r) reaches its maximum 
(inset), e, In the absence of experimental noise, for each triplet of birds, i,j, k we must have, Tik + T^j = Tij\ if i turns 20ms 
before k, and k turns 15ms before j, then i turns 35ms before j (Time Ordering Relation - TOR). Due to noise TOR will not 
hold strictly, but we still want it to be correct on average for Tij to make biological sense. We consider all triplets of birds and 
plot + Tkj vs. Tij. The data fall on the identity line with relatively small spread, confirming the temporal consistency of 
the turning delays. 



that the equations commonly adopted to describe collec- 
tive motion fail to explain. We introduce a new theory 
mathematically identical to that describing superfluid- 
ity in liquid helium. The superfluid theory not only 
explains the observed sound-like propagation, but also 
predicts that the speed of propagation must be larger in 
more ordered flocks, a prediction confirmed by the exper- 
imental data. This result indicates that the efficiency of 
information transport during collective decision-making 
is quantitatively linked to the degree of behavioural po- 
larization in the group. 

We study natural flocks of starlings performing aerial 
display at dusk [T2 [T3]. By using a 3-cameras setup 
we reconstruct the full 3d dynamical trajectory of each 
bird in the flock. We have negligible time fragmentation: 
90% of the reconstructed trajectories last more than 90% 
of the duration of the studied event. Cameras shoot 
at 170Hz (Methods). From our pool of data we select 
12 flocks, each one performing a collective turn. Turns 
typically last a few seconds. In Fig. [T]we present samples 
of the reconstructed trajectories. When bird i makes a 
turn, the modulus of its radial acceleration, Oj(i), has a 
maximum (Fig. [l]-d) . We exploit this simple kinematic 
fact to organize all birds in the flock according to their 
temporal relationships. For each pair of birds, i and j, 



we calculate their mutual turning delay, , namely the 
amount of time by which bird j turns before (ry > 0) 
or after (r^ < 0) bird i (Fig. Illd,e) . By using the delays 
t^ we then rank all birds in the flock according to their 
turning order, that is we find who is the first to turn, 
who is second, and so on. In this ranking, each bird i is 
labelled by its rank, r i; and by its absolute turning time, 
ti, which is defined as the delay with respect to the top 
bird in the rank, i.e. the first to turn (see Methods and 
Appendix [F] for details) . To represent the ranking in a 
compelling way, we plot the rank of each bird as a 
function of its absolute turning time ti, thus obtaining 
the ranking curve, r(t) reported in Fig. [2ja. 

By using the ranking we can find out whether the turn 
has extended or localized spatial origin. In the extended 
case, a large number of birds start to turn all at the same 
time, hence we would expect to find many birds packed 
in a very short time lag at the top of the ranking, i.e. 
a negative second derivative of r(t) for small t. What 
we see from the data is the opposite: the ranking curve 
r(t) is convex for early times, meaning that the turn is 
started by very few birds (Fig. |2]-a). Moreover, we find 
that the first birds to turn (say, the top 5 in the rank) 
are physically close to each other. More precisely, their 
average mutual distance D does not scale with the size 
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FIG. 2. Propagation of the turn across the flock, a, The rank r of each bird in the flock, i.e. its order in the turning 
sequence, is plotted vs its absolute turning delay t, i.e. the delay with respect to the top bird in the rank (the first to turn). 
The convex toe of the curve for early times indicates that few birds take initially the decision to turn, b, The average mutual 
distance D between the top 5 birds in the rank does not increase with the linear size of the flock, L, hence indicating that the 
first birds to turn are actually close to each other in space. The result does not change if we use a different number of top 
birds. Inset: the actual position of the top 5 birds (red) within a real flock, c, The distance x traveled by the information in 
a time t is proportional to the radius of the sphere containing the first r(t) birds in the rank, namely x(t) = [r(t)/ p] 1 ^ 3 . The 
linear regime of x(t) allows us to define a 'sound' speed of propagation, c s , as the slope of x(t) for early-intermediate times. 
The speed c s varies significantly from flock to flock (see also Table EJ. d, The intensity of the peak of the radial acceleration, 
a max , (solid symbols) decreases very weakly in passing from the first to the last turning birds. In the inset, we plot af ax vs the 
rank r; for each bird. This slow decay indicates that the information propagates through the flock with negligible attenuation. 



L of the flock (Fig. [2]-b) . This result indicates that the 
number of birds initiating the turn is not proportional 
to the volume {D ~ L 3 ), nor to the surface (D ~ L 2 ) of 
the flock. Hence, the decision to turn has a spatially lo- 
calized origin and it then travels across the flock through 
a transfer of information from bird to bird. 

To understand how effective is this transport of infor- 
mation we need to calculate how much distance x the 
information travels in a time t, i.e. we need the dis- 
persion law. We are in three dimensions and the turn 
has a localized origin, hence x(t) is equal to the radius 
of the sphere containing the first r(t) birds in the rank, 
namely x(t) — [r(t)/ p] 1 ^ 3 , where p is the density of the 
flock (Fig. [2j-c) . The most striking feature of the prop- 
agation curve x(t) is that there is a clear linear regime 
for early-intermediate times (before border effects kick 
in). The distance traveled by the information grows lin- 
early with time, x(t) = c s t, just as a sound wave. The 
parameter c s is the speed of propagation of directional 
information, which is in the range 10 — 20 meters per sec- 
ond (Table I). The second important result is that the 
information to turn propagates across the flock with neg- 
ligible attenuation (Fig. [2]-d). This is nontrivial: flocks 
are large, the information to turn dynamically reaches 
all birds through a lot of intermediate passages, so that 



a substantial level of damping could be expected. Yet it 
is not so. This phenomenon too is reminiscent of sound 
propagation. 

The speed of propagation of the information, c s , varies 
significantly from flock to flock (Fig. [2]-c and Table EJ. It 
therefore seems that some flocks are more efficient than 
others to transport information. Why is that? Thinking 
about sound, one may naively expect the speed to de- 
pend on density. However, the variability of c s does not 
disappear by rescaling it with the flock's density. In fact, 
even though a linear dispersion law is suggestive of sound 
propagation, we should not forget that what propagates 
during the turn are fluctuations of orientation, not of 
density. We shall see later that the variability of c s has 
an entirely different explanation. 

Linear propagation and low damping of the signal are 
key factors in achieving an efficient collective decision 
of the group. Both sub-linear propagation and atten- 
uation would result into a physical spread of the flock, 
and eventually into total disruption. Do current theories 
of collective motion account for such an efficient trans- 
port of directional information? Virtually all theoretical 
descriptions are based on alignment dynamics: each in- 
dividual tends to keep its direction of motion as close as 



4 



possible to that of its neighbours 



i(t + 1) = Vi(t) + Jj^v^t) , 

jei 



(1) 



where the vector Vj is the velocity of bird i and 
the sum extends over all neighbours j of i (be they 
metric or topological [T2J, [23]). We have disregarded 
noise/temperature, which is inessential for what follows; 
we will just assume that the alignment strength J is 
large, so that we are in the deeply ordered phase (as 
natural flocks are |16j). In continuous time equation ([I]) 
is equivalent to a zero-temperature Langevin equation, 



dt 



dH 
dvi 



H 



-J^Vi-Vj 



(2) 



According to these relations each bird updates its veloc- 
ity following a social force, F* = —dH/dvi, produced by 
its neighbours. The Hamiltonian H in ^ is the same as 
that of a ferromagnetic system, where the birds velocities 
Vj play the role of magnetic spins [24j . 

To simplify the algebra we exploit the fact that the 
trajectories of birds during a turn lie approximately 
on a plane (Fig. [T]-b,c). This allows us to use a two- 
dimensional velocity, = (v'f,vf) — ve lipi , where the 
phase ifi is the angle between the direction of motion of 
i and that of the flock (we make the standard assump- 
tion that v is constant) . In the highly ordered phase the 
velocities differ little from the collective one, so that 
PiCl We can thus expand H in eq.Q [23] , 



H = g ~ ^J') 2 = \ a2j 



^ [Vp>{x,t)f , (3) 



where a is the average nearest neighbours distance and 
a term v 2 has been reabsorbed into J. The Langevin 
equation associated to Hamiltonian ^ is, 



dtp 
~dt 



SH 

Sip 



a 2 JV 2 ip 



(4) 



Relation Q is a diffusion equation for the phase ip, and 
it has dispersion law to — ik 2 . This result has two conse- 
quences, both in sharp contrast with the empirical data: 
i) information travels much slower than linearly, x ~ vt, 
at variance with the linear propagation we find in turning 
flocks; ii) the frequency is imaginary, meaning that this 
is a non-propagating mode. Transfer of information gets 
damped exponentially in space and time, again in stark 
disagreement with the brisk, undamped propagation we 
observe in flocks. 

The standard theory has two problems. First, it seems 
to be missing some conservation law. Hamiltonian ([2| is 
invariant under a global gauge symmetry, namely the 
uniform rotation of the velocities v, (ifi — > ifi + Sip). 
This symmetry encodes the fact that all directions of 
flight are equivalent for a flock. Through Noether's theo- 
rem, a symmetry implies in general a conservation law, of 



which, however, there is no trace in the standard theory. 
A hidden conservation may heavily affect the dispersion 
law, because a conserved quantity cannot be relaxed lo- 
cally, but it must be transported across the system. Sec- 
ond, equation Q completely neglects behavioural iner- 
tia, as the social force, F s = aJW 2 ip, controls directly 
tp, rather than ip. This is odd: imagine that the inter- 
action with the neighbours requires bird i to perform a 
U-turn in one time step. This behaviour is allowed by 
the standard theory, although it is clearly unreasonable. 

To address these problems we follow Landau's ap- 
proach [31] , namely we identify a suitable order parame- 
ter and write the simplest Hamiltonian compatible with 
the symmetries and the constraints of the system. The 
phase ip is the obvious order parameter. The gauge sym- 
metry implies that H can only be a function of V</?, not 
of ip itself, as expressed by ([3]). On the other hand, the 
inertial constraint requires a kinetic term, s 2 /2\, where 
s z is the canonical momentum conjugated to ip, and \ is 
the generalized moment of inertia. We therefore propose 
the novel Hamiltonian, 
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(5) 



where p s = a 2 J, is the rescaled alignment coupling con- 
stant, or stiffness [30 . The momentum s z is defined 
as the local generator of the rotations parametrized by 
the phase tp., so that (s z ,ip) are generalized action-angle 
canonical variables. It can be shown that s z is essentially 
the inverse radius of curvature of the trajectory, whereas 
X is the behavioural resistance of a bird to change its in- 
stantaneous radius of curvature when a social force is 
exerted by its neighbours (see Appendix [A] and [B]) . The 
canonical equations of motion generated by ([5]) are, 



dip 
~dt 



S_H_ 

5s z 



X 



ds z 
dt 



5_H_ 

Sip 



p s V ip 



(6) 



The crucial consequence of the gauge symmetry is that 
the r.h.s of the second equation of motion is in fact a 
gradient, V 2 ip = V • V<y3, so that we can rewrite this 
relation as a continuity equation for s z (x,t), 



ds z 
dt 



Vj 2 







(7) 



with current j z (x,t) = p s Vp{x,t). We therefore find a 
conservation law. Imagine that a strong misalignment 
among a subgroup of birds forms in a certain position 
of the flock. This causes a local excess of curvature, 
and thus an excitation of the field s z {x,t). Conserva- 
tion law (|7|) states that such excitation cannot be locally 
dissipated out, but it must be transported away. This 
mechanism, which is the essence of the new theory, gives 
rise to an undamped sound-like mode. Indeed, by taking 
the second derivative with respect to time in ^ we get, 



dt 2 



Ps/X 



(8) 
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FIG. 3. Superfluid prediction, a, The superfluid theory predicts that the rescaled speed of propagation of the turn, c a /a, 
must be a linear function of l/yl — where $ is the polarization. The prediction is verified by the empirical data (P- value: 
P = 3.1 x 1CP 4 ; correlation coefficient: R 2 = 0.74). Each point is a different turning flock. Error bars on c s are obtained 
from its variability under changing the linear fitting regime of x(t). c 3 /a has the dimensions of sec -1 . The slope of this line is 
equal to \ftjx ~ equation) 10 1, b, Polarization as a function of time in three different flocks. The value of $ reported in panel 



a corresponds to the time average over the entire duration of the turn. 



Relation ^ is D'Alembert's equation, describing waves 
propagating with speed c s and no damping. Its disper- 
sion relation is linear, uj = c s k, which in terms of the 
distance x traveled by the information in a time t, reads, 
x = c s t. This is precisely the linear and undamped prop- 
agation law that we find in turning flocks. 

Some flocking theories have established a correspon- 
dence between flocks and magnets, with birds veloci- 
ties v, playing the role of spins trying to align to each 
other |19l 1201 124j . However, up to now spins were vir- 
tual, as they did not obey proper Poisson rules. Within 
the present description, things change. The fact that s z 
generates the symmetry parametrized by the phase tp is 
expressed by Poisson relation, {v, s z } = dv/dip = iv, 
which in components reads, 



ii) spontaneous symmetry breaking, i.e. emergence of 



{v x ,s z } 



dv x 
dp 



{Vy,S z } 



dp 



(9) 



If we interpret v x and v v as the x, y components of the 
spin, equations ^ show that s z is a true spin, namely 
the generator of the rotation in the space of the order 
parameter v. This is the most general and fundamental 
definition of spin [35]. Accordingly, ] z = pgV(p is the 
spin current and p s the spin stiffness |34l 135] . 

The theory of collective motion that we have intro- 
duced above is exactly the same as that describing su- 
perfluid liquid helium (He-II). It has been demonstrated 
long ago by Matsubara and Matsuda [32j [33] that the 
lattice-gas model for Bose condensation in He-II is math- 
ematically equivalent to the planar ferromagnetic model 
defined by equation (j5J). This is not just a bizarre 
coincidence. Superfluidity is nestled into the identical 
mathematical structure of these apparently very differ- 
ent systems. The keystones of superfluidity are: i) ex- 
istence of a gauge symmetry (arbitrariness of the quan- 
tum phase/arbitrariness of the flock direction of motion); 



a nonzero order parameter, xj) 



j e lv (nonzero Bose 



wave function/nonzero flock's velocity); iii) coupling of 
the phase ip to the generator of the gauge symmetry 
s z (Bose particle density/spin density). Irrespective of 
the physical and biological details, these three elements 
alone generate the dissipationless propagating mode de- 
scribed by equation (JsJ) , that is superfluidity [3~4Tf3"T] . In 
liquid helium, j 2 is the current of the superfluid compo- 
nent and the propagating mode is called 'second-sound' 
[34] [38]. In flocks, j z transports spin, that is curvature, 
giving rise to the collective turn. 

The superfluid theory not only provides an explana- 
tion for the linear and undamped propagation of infor- 
mation in natural flocks, but it also makes a prediction 
about the dependence of the speed of propagation on the 
experimentally accessible quantities, thus making sense 
of the otherwise unexplained variability of c s from flock 
to flock. We recall that c 2 s = p s /x> where the stiffness 
is p s = a 2 J. The alignment strength can be written as 
J = e/(l — <!>), where <& is the polarization, namely the 
degree of alignment in the flock, $ = ||(1/JV) J2i v i/ u ill> 
and e is an energy constant, setting the scale of the align- 
ment interaction |2"4"1. We therefore obtain, 



(10) 



Equation ( 10 1 predicts that the speed of propagation of 



the turn across a given flock must be larger the larger 
the degree of alignment <!> in that flock. We report c s /a 
vs. 1/Vl — $ for all our flocks in Fig. [3] Data show a 
clear linear dependence, exactly as predicted by equation 
( |10p . We remark that the square root behaviour repro- 
duced by the data is nontrivial: the polarization is a 
dimensionless quantity, hence the functional dependence 
of c s on $ cannot be worked out by mere dimensional 
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analysis. Interestingly, equation ( 10 ) also implies that 



the only empirically observable quantity is the ratio be- 
tween the scale of the alignment interaction, e, and the 
turning inertia x, so that different species (or artificial 
entities) may have very different values of e and %, but 
still have a comparable information transport efficiency. 

Superfluid propagation of information in collective 
motion has never been discovered before, not even in 
models more realistic than 0, as those studied in [251 - 
[25] , In fact, as we have seen, superfluidity stems from 
associating inertial terms to the fundamental conserva- 
tion law generated by the gauge symmetry. It would be 
interesting to reconsider flocking models with inertia in 
the light of the present results. Moreover, we note that, 
although the dynamical equation Q is different from 
Q, the static properties of the field (p(x,t), and in par- 
ticular its equal time correlations are the same as 
those described by Hamiltonian (|3). This result is due 
to the separation between coordinates and momenta. 

The link between speed of propagation of the infor- 
mation, c s , and behavioural polarization, $, is not an 
evolutionary trait, but the mathematical consequence of 
the gauge symmetry. However, the specific level of polar- 
ization of a flock is not fixed by math, nor by symmetry, 
but by adaptive factors. In many social species polar- 
ization is very large [H [5j [3 [16] . Global order is indeed 
the most conspicuous trait of collective behaviour. How- 
ever, were the only concern of a bird not to bump into 
its neighbours, such a large polarization would be diffi- 
cult to justify. Flocks for which we have data are rather 
diluted systems, with packing fraction lower than 0.01 
[T5] . Yet these same flocks are very stiff, precisely in the 
superfluid sense: alignment strength, p s ~ J, is large 
and polarization is close to 1. Why is that? In collective 
decision-making swift transport of information is bene- 
ficial to the group. In the case of turns this is obvious: 
during the turn the wavefront divides approximately the 
flock into two groups of birds with different directions 
of motion. Such misalignment causes a spatial spread of 
the flock, with overall loss of cohesion. The slower the 
speed c s of the wavefront, the more severe this loss. It is 
reasonable to believe that this is a general mechanism. 
Every collective decision drives the group through a mo- 
mentary lapse of cohesion, due to the transient coexis- 
tence of different behavioural states. The superfluid link 
between high behavioural polarization and fast propaga- 
tion of the information suggests that keeping this lapse 
to a minimum, therefore achieving a swift and robust col- 
lective decision, may be the adaptive drive for the high 
degree of order observed in living groups. 
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Appendix A: The spin 

Let us define external space the space of the birds co- 
ordinates and internal space (or target space) the space 
of the order parameter, namely the planar velocity. In- 
deed, v(x, t) is a map between the external space E 3 x K 
and the internal space SO (2) (the circle). It is essential 
to understand that the phase ip parametrizes rotations 
in the internal space of the velocity and it must not be 
confused with the standard angle 6 of 2d polar coordi- 
nates, which parametrizes rotations in the external space 
of positions. The easiest way to understand this is the 
following: a bird flying straight (with respect to an ar- 
bitrary fixed reference frame - that of our cameras, for 
example) has <fi = 0, but 6 ^ 0. 

It is interesting to note that rotations parametrized by 
these two angles, tp and 8, correspond to two very dif- 
ferent types of collective turns. A rotation parametrized 
by ip corresponds to an equal radius turn, i.e. a turn 
in which all birds have (approximately) the same radius 
of curvature and where trajectories cross. This is ex- 
actly the way adopted by flocks to turn (see Fig. [I]), 
and it has a clear biological motivation: it keeps the 
speed v (approximately) constant throughout the flock. 
Equal radius turning was first experimentally discovered 
in [32] (see also [T3]). On the other hand, a rotation 
parametrized by 9 corresponds to a parallel path turn, 
typical of rigid bodies. This kind of turn implies differ- 
ent radii of curvature for different points, and therefore 
it requires an increase of speed of birds on the outer side 
of the turn, which is clearly biologically unreasonable. 

The generator of the external ^-rotations is the stan- 
dard angular momentum, l z , whereas the generator of 
the internal (^-rotations is the spin, s z , which is con- 
served by the continuity equation. In order to have some 
intuition about the physical nature of s z we must con- 
nect external to internal spaces. This connection is es- 
tablished by the kinematic equation, 



v e 



(Al) 



expressing the simple fact that birds are not anchored 
to a lattice, but they follow their velocity vectors. If we 
cons ider the speed v approximately constant, equation 
( Al I implies, 



ip = v/R 



(A2) 



where R is the instantaneous radius of curvature. Using 
(A2| together with we get, 



v ■ s z 



(A3) 



and dividing by the first equation in ^ we finally obtain, 

(A4) 



s z R 
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This equation connects the rather abstract internal spin, 
s 2 , to a very clear kinematic quantity, the radius of cur- 



be rewritten as, 



vature R. By integrating (A4), we have 



l/R 



(A5) 



where k — l/R is, by definition, the curvature of the 
trajectory. Therefore, once the connection with the ex- 
ternal space is performed, the spin turns out to be essen- 
tially the curvature. This is why a bird flying straight 
(R = oo, K ~ 0) has s z ~ 0, while it has nonzero stan- 
dard angular momentum l z . A change (in time) of the 
spin s\ of bird i, due to the social force exerted by the 
neighbours of i, corresponds to a change (in time) of its 
instantaneous radius of curvature, Ri, and curvature, «j. 
Hence, what actually propagates across the flock during 
the turn is a fluctuation (in space and time) of the cur- 
vature field, K,(x,t). Before the turn, the flock is flying 
almost straight, R 3> 1, k ~ 0, s z ~ 0. Then the turns 
sparks in some part of the flock, causing an increase of 
the curvature k, i.e. an increase of s z . This change 
sweeps through space and time until the whole flock has 
turned. Finally, after the turn, the flock relaxes back to 
R ^> 1, k ~ 0, s z ~ 0. Mathematically, this propagation 
of the curvature, i.e. of s z , derives from the canonical 
equations ([6]): by taking the second derivative with re- 
spect to time one obtains a D'Alembert wave equation 
for s z {x,t) identical to that obeyed by <p(x,t), eq.Q. 

Conservation law ([7]) states that the spin-curvature, 
s z (x,t) <~ n(x,t), obeys a continuity equation. As we 
have seen, this conservation law is crucial to determine 
sound-like propagation. Continuity means that the tra- 
jectory curvature in a given volume cell of the system 
cannot change unless it is transported into, or out of, 
that cell by a spin current, j 2 = p s Vip. We can reformu- 
late this by saying that, if a certain excess of curvature, 
namely a strong misalignment among a certain group of 
individuals, forms in a given point of the system, it can- 
not be simply dissipated out. Rather, such excitation 
creates a social force that makes the neighbours turn, 
and their neighbours too, and so on, so that the excess 
curvature is transported away, instead of dissipated. 

Finally, note that the spin s z is not the z component 
of the velocity. Also note that the rotation generated by 
s z is the very transformation under which Hamiltonian 
([5]) is symmetric. 



Appendix B: The generalized moment of inertia 

As we have seen, the spin s z is not the standard angu- 
lar momentum l z . Accordingly, the generalized moment 
of inertia x, is not the standard moment of inertia, which 
in the case of circular motion is, I = mR 2 , where m is 
the mass. So, what is the physical and biological mean- 
ing of x? From the canonical point, the answer is clear: 
X is the inertia to changing ip. Indeed, equation ([8| can 



X 



aF s 



(Bl) 



where F s — aJV 2 tp is the social force exerted by the 
neighbours. Hence, the generalized moment of inertia x 
is defined as the ratio between the social force (the cause) 
and the change of angular velocity (the effect). This is 
the standard definition of inertia: the ratio between force 
(cause) and acceleration (effect). However, in this case 
F s is a generalized (social) force, and aF s is a generalized 
torque, hence x is not the standard moment of inertia. 

To better grasp the biological meaning of x we must, 
once again, bridge the gap between internal and exter- 
nal space. By using equations (Al), (A2) and (Bl| we 



obtain several cause-effect relations shading light on the 
physical and biological meaning of x- The first relation 
connects the social force to the change of radius R, or 
cquivalently to the change of curvature ft, 



X 



v 



aF s 
R 



aF s 
k 



(B2) 



Hence, x is the resistance of a bird to change its instan- 
taneous radius of curvature R (the effect), when a given 
social force F s (the cause) is exerted. 

Another interesting relation can be obtained in terms 
of the banking angle 7. A banked turn [39l [40] is the 
typical way birds (and planes) change their heading. It 
consists in a gentle roll, so to form an angle 7 between 
the axis of the wings and the horizontal plane. In this 
way, part of the total lift goes into a centripetal force, 
F c = mg"f, making the bird turn (m is the mass, g the 
gravitational acceleration and 7 <§C 1). From (B2), it is 
straightforward to prove that, 



aF s 

7 



(B3) 



According to ( B3 ) , the generalized moment of inertia x 



is the resistance of a bird to change its banking angle 7. 
Here, 7 is the effect of the social force, F s , and x se ^ s 
the ratio between cause and effect. Notice that equations 



( B2 ) and ( B3 ) arc clearly non-canonical definitions of the 



inertia x, because at the denominator they both have a 
first order derivative in time, rather than a second order 



one, as in the canonical equation (Bl ) 



Let us emphasize once again that x 1S n °t the stan- 
dard, mechanical moment of inertia, I = mR 2 , but 
rather a social, or sensorimotor, resistance of the bird 
to change R or 7. It is not possible to write an ab inicio 
expression for x m terms of primary mechanical quanti- 
ties, like mass, radius, etc. To understand this fact, let 
us imagine that at some point the neighbours of bird i 
all sharpen their banking angle 7, thus creating a strong 
social force, F s = aJV 2 ip, acting on i. What we call a 
social 'force' is in fact a shortcut to describe a very com- 
plex sensorimotor process: a nonzero V 2 <p means that 
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some of the neighbours of i are now about to crash into 
i. This is most likely perceived by i, which decides to 
change its own 7 and make it equal to that of the neigh- 
bours, thus avoiding the crash. However, the degree by 
which i will react to the imminent crash, or conversely 
the resistance to this reaction (which is x), is the result 
of a very complex trade-off. Let us analyze this trade-off 
by pretending to be i. 

On one side, there is the price of the crash. How 
imminent is it? This will depend on both the nearest 
neighbour distance and on the mutual velocity. How 
bad would that be? Perhaps, I can ignore my neigh- 
bours, and just change them, without any real crash. 
How confident I am into my capability to change 7? If I 
am very agile, I can wait a bit longer before changing my 
7. On the other side, there is the price to changing 7. 
How much will it cost me energetically? By increasing 7 
I will increase the drag against air, otherwise I fall down. 
But to do this I must increase the power, which is very 
costly. Can I manage? 

The generalized moment of inertia \ is the very final 
product of this process. Clearly, we cannot know a priori 
its value. But we can define it, and measure it. This is 
exactly what we have done by measuring the speed of 
propagation of the turn across the flock. In this sense, 
the situation is the same as in real magnetic systems 
[34] : the parameter x 1S the magnetic susceptibility to 
an external field, which cannot be simply expressed as 
a function of the microscopic parameters of the theory, 
but it can be experimentally measured. 



Appendix C: Dimensional analysis 

In the supcrfluid theory we have an Hamiltonian that 
is the sum of two parts. Hence, we have to be care- 
ful with physical dimensions. The phase is of course a 
pure number, [<p] = [1], whereas the alignment coupling 
constant has the dimensions of an energy, [J] = [e]. In 
this way the social force has the dimensions of a true 



force, [F s ] — [aJV 2 <p] = 



_1 ] and the spin has the 



dimensions of an angular momentum, [s z ] = [e ■ t], i.e. 
of an action. Accordingly, x nas the dimensions of a 
true moment of inertia, [x] — [e • t 2 ]. Notice that the 
term appearing in equation ^ is a 2 JV 2 ip = aF s , which 
is dimensionally a torque. Hence, the derivative of an 
angular momentum is a torque, as it should. 

By definition, the polarization is a pure number, [$] = 
[1]. This is why in the relation linking alignment cou- 
pling constant to polarization, J = e/(l — $), we need 
a dimensional constant with the dimensions of an en- 
ergy, [e] = [e]. As we have said, e sets the scale of the 
alignment interaction. Roughly speaking, e (as well as 
x) is what distinguishes species A from species B, or 
bird flocks from aircraft formations. Finally, with to the 
above physical dimensions, the speed of propagation of 
information across the flock, c s , is measured in meters 
per second, as it should. 



Appendix D: General off-plane case 

Our initial assumption that the birds' velocities lie 
on a plane during the turn, namely that the turn has 
very small torsion, although experimentally satisfied (see 
Fig. lb, c), it is not at all a necessary condition for our 
mathematical description. The most general formulation 
of our result holds even with a truly 3d order parameter 

v 2 G2|. 

If we assume that the mean direction of motion of the 
flock points in the x direction, then there will be full 
3d fluctuations of the individual velocities Vj around the 
overall direction of motion of the flock, generating small 
components of Vj along the two orthogonal axes, z and 
y. Therefore, we must define two phases, ip z and ip y and 
the Hamiltonian can be spin-wave expanded in terms of 
these two fields. The phase (p z parametrizes rotations 
of Vj around the z axis (as in the planar - zero torsion 
case), whereas tp y parametrizes rotations of Vj around 
the y axis. In this fully 3c? case the Hamiltonian is given 
by [51, 



H = 



d 3 x 1 



(V^ 2 ) 2 + (V^) 2 



1 

2x 



si + s y ] 



where p s — a 2 J is, as usual, the stiffness. The equations 
of motion are, 



dtp a 


dH 


dt 


ds a 


ds a 


dH 


dt 





X 

= a 2 JV 2 tp a = V • j Q 



(D2) 
(D3) 



with a = y, z, giving rise to two D Alembert equations, 



d 2 (f a 
dt 2 



c 2 s VVa = 



Pslx 



(D4) 



In the full 3d case we therefore obtain two, rather than 
one, propagating dissipationlcss modes, along the trans- 
verse directions y and z. This is just a manifestation of 
Goldstone's theorem [4"4"j . 

These equations are exactly the same as for model 
G in the Halperin-Hohenberg classification of dynami- 
cal universality classes [3H [35] . Model G does not de- 
scribe supcrfluid liquid helium, but an isotropic Heisen- 
berg antiferromagnet with staggered magnetization as a 
non-conserved order parameter, and total magnetization 
as a constant of motion. An experimental realization of 
a 3d isotropic antiferromagnet is RbMnF3, a compound 
whose dynamics is characterized by the transverse spin- 
wave modes ( D4 1 . Notice that also in this system there 
is superfluid transport. As discussed in the main text, 
superfluidity is not restricted to liquid helium II, but it is 
rather built into the mathematical details of the theory. 
In particular, it is the product of symmetry and conser- 
vation laws. In the full 3d case described here (model G) 
these ingredients give rise to superfluid transport exactly 
as in the planar (He-II) case. 
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To write (Dl ) and (D4| we have assumed that the two 



excitations <p z and <p y are equally likely, so that the only 
symmetry breaking direction is that of motion. In fact, 
recent studies on individual diffusion in starling flocks 
show that gravity is another symmetry breaking direc- 
tion, heavily suppressing fluctuations along the vertical 
plane |45j . If we identify z with the axis of gravity, this 
suppression would imply that rotations of the velocity 
around the y axis are suppressed, and therefore that (p y 
is less relevant a degree of freedom than <p g . This sup- 
pression induced by gravity is likely the cause of the 
planar-like turns we observe in flocks and it thus justi- 
fies the adoption of the simpler planar description of the 
main text. 



With small dissipation there is linear propagation of the 
information and the time scale of the exponential decay 
is set by r = 2\/r]. From (E6| we get, 



> Vx/Ps L = L/c s 



(E7) 



Therefore, small dissipation implies that the damping 
time constant is larger than the time the information em- 
ploys to travel across the flock. In other words, the signal 
is effectively very weakly damped across the length scale 
of interest. We conclude that even when a small dissi- 
pation is present, propagation of information is qualita- 
tively the same as that described by the zero dissipation 
theory. 



Appendix E: Dissipation 

The diffusive equation of motion Q derived from the 
standard theory can be obtained as the over-damped 
limit of the new theory, once we introduce a dissipative 
term proportional to ip in the equation of motion, 

with p s — a 2 J and where 77 is a generalized friction co- 
efficient. From this we get, 



XOJ — ITjLO — p s k 







(E2) 



In the limit r\ 3> X we simply get the diffusive result, 
uj = i(p s /rj)k 2 . In general, however, we obtain, 



ui 



1 — ± 
2X 



c s k\J\ k 2 /k 2 , 



(E3) 



where, as usual, the propagation speed is c s = \fp~7Jx 
and, 



fco 



V 



(E4) 



If we define the dissipation time scale, r = 2x/rj, and 
the zero-dissipation frequency, ujq = c s fc, we can rewrite 
the dispersion law as, 



uj = i Jt ± ujq yjl — k^/k 2 



(E5) 



With zero dissipation, we get fco = 0, r = 00 and uj = ujq. 
For r\ 7^ 0, on the other hand, we have two regimes, ac- 
cording to the value of the friction coefficient and of the 
wave number k. For k > fco we have attenuated propa- 
gating waves, as the frequency has both a real and an 
imaginary part. For k < fco we have evanescent waves: 
the frequency is purely imaginary, there is no propaga- 
tion, but pure exponential decay. 

The smallest value of fc in the system is fc m i n ~ 1/L. 
where L is the linear size of the flock. Hence, small 
dissipation is defined by the relation, 



small dissipation 



(E6) 



Appendix F: Mutual delay vs reaction time 

One may think that the mutual delay between two 
birds, Tij, is the same as the reaction time, tr, namely 
the time between the stimulus provided by j and the 
consequential action of i. However, this is not the case. 

Let us assume j is the first bird to turn, and that i is 
second. By definition of reaction time, i begins its turn 
tr seconds after j. However, we do not define as the 
difference between the starting instants of the two turns: 
there is no practical and robust way to do that, because 
each birds turns smoothly, so that it is impossible to 
define the 'start' of the turn. To compute t, 7 we use 
the entire trajectory of both birds, by finding the time 
shift that maximally overlaps the accelerations of i and 
j (see Fig.l). If the function ai(t) were exactly the same 
as the function a,j(t — Tij), then we would have r^- = tr. 
This, however, is never the case. First of all there is 
noise, so that the two curves, <Xj(i) and Oj(t), are only 
approximately shifted with respect to each other. But 
more importantly, the second bird j can try to 'catch up' 
during the turn, so that the delay at the end of the turn 
is shorter than the delay at the beginning of the turn, 
which is the reaction time. In this case, the delay T.y 
would be a value intermediate between those two times, 
hence giving a value smaller than tr. The opposite can 
happen too: bird i could in fact lose ground during the 
turn, so that the delay at the end of the turn is longer 
than the reaction time, and r,j is larger than tr. 

There is, however, a robust way to extract a timescale 
from our data, and to compare it to standard reaction 
times in birds. The speed of propagation, c s , is expressed 
in meters per second, so that c s /a is the inverse of a 
time constant: it is essentially the time the information 
needs to travel the nearest neighbour distance, a. From 
Fig. 3 we see that a/c s ranges between 25ms up to 100ms, 
with an average around 50ms. Again, this is not the 
reaction time, as it depends on the information transfer 
mechanisms. However, it is reassuring to see that 50ms 
is in the physiological range of reaction times for birds 
in general and for starlings in particular 43 . 
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FIG. 4. Check of time ordering relation (TOR). We report the TOR test for several of our flocks and for one non-turning 
flock (lowest-right panel) . Temporal consistency requires that Tij ~ Tj& + rtj , so to have the data scattered along the identity 
line with clear correlation. In the case of the non-turning flock, on the other hand, the delays are just random numbers, so no 
temporal consistency is found. 



METHODS 



Experiments. European starlings (Sturnus vulgaris) 
spend the winter in Rome, where they populate several 
roosting sites. Data were collected at the site of Piaz- 
za dei Cinquecento, between November 2010 and De- 
cember 2012. To acquire the video sequences we use 
the trifocal stereometric setup described in [14 . Two 
cameras separated by a baseline distance D\i = 25m 
are the stereometric pair. A third camera, placed at 
a shorter distance D23 — 2.5m from the first one, al- 
lows us to exploit the trifocal constraint for solving the 
stereo correspondence (matching) problem [14] . We em- 
ploy three high-speed cameras IDT-M5 with monochro- 
matic CMOS sensor with resolution 2288 x 1728 pixels, 
shooting at 170hz. Cameras store images on off-board 
memory using the Camera Link protocol. Lenses used 
are Schneider Xenoplan 28mm f /2.O. Typical exposure 
parameters are: aperture between f/2.8 and f/8; ex- 
posure time between 700 and 3500 ms. Intrinsic camera 
parameters are calibrated in the lab using a set of images 
of a known planar target. The recorded events have a 
time duration between 1.8 and 12.9 seconds. The data- 
set consists of 12 distinct flocking events, each one in- 
cluding one collective turn. If on the recorded sequence 
the flock performs more than one turn, the time lag is 
chopped and different turns are studied as independent 
events. 



Tracking. Detection of individual birds on the im- 
ages is carried out by using the same method as in 
[14]. To assign stereoscopic links, i.e. to match birds 
across the three images, we use global optimization us- 
ing a cost function based on the trifocal constraint [46] , 
To perform temporal linking we first use a roto-scale- 
translation |47j to predict the position of each bird in 
the next frame. We then link birds from one frame to 
the next one in a redundant way, i.e. when in doubt we 
use multifurcation. Percolating the full set of temporal 
links though the entire temporal sequence, we construct 
all the possible 2d paths in the image space of each one 
of the three cameras. The three sets of 2d paths are then 
matched via a global assignment, based on a cost func- 
tion which measures the number of stereoscopic links be- 
tween each triplet of 2c? paths. Global optimization not 
only matches the right 2c? paths, producing the correct 
3c? trajectories, but also eliminates the unphysical paths. 
To avoid exponential explosion of the number of paths, 
the temporal sequence is recursively divided into shorter 
time intervals over which the complexity of global opti- 
mization can be handled. All global optimizations are 
performed using linear programming |48j . Full details of 
the tracking algorithm will be reported elsewhere. 

Filtering. Filtering of the time-discrete trajectories 
is necessary for two reasons: i) to reduce experimen- 
tal noise; ii) to eliminate wing-flapping, whose frequency 
for starlings is Wfl ap = lOhz. By sampling at 170hz we 
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TABLE I. Polarization and speed. N is the number 
of birds in the flock. The polarization is defined as, <& = 
||(1/AT) Vi/^ijj. The values of $ reported here are on av- 
erage smaller than those reported in previous investigations 
[121 113] . The reason for this is that previous data were ob- 
tained with cameras sampling at 10Hz, whereas the present 
data are obtained at 170Hz. At this sampling frequency ex- 
perimental noise and wing flapping reduce the polarization. 
This reduction, however, affects equally all flocks by simply 
rescaling (1 — hence it does not change the correlation 
m Fig(3] The speed of propagation of the information, c s , is 
found by fitting the linear regime of the propagation curve, 
x(t). We note that by fitting the ranking curve to a power 
law for early-intermediate times, r(t) — t a , we find on average 
a = 3.2, thus justifying the statement that x ~ r 1 ^ 3 ~ t 1 ' 07 
is a linear function. 



are fully exposed to the trajectories zig-zag (see inset in 
Figjl^,), which would completely dominate acceleration. 
To cut this high frequency mode we a use 2nd order low- 
pass digital Butterworth filter on the velocities, typically 
with a cutoff frequency Wfl ap /30. Accelerations obtained 
in this way (see FigjTJi) capture the low frequency direc- 
tional change corresponding to the turn. Our final re- 
sults are robust against changes of the cutoff frequency. 

Turning delay. We define the turning delay Ty of 
bird i with respect to bird j as the time by which we 
have to shift the radial acceleration a jit) with respect 
to diif) to maximally overlap them. More precisely, we 
define the following normalized correlation (or overlap) 
function, 



Given this definition, the time shift Tij corresponds to 
the value of r where Gy(r) reaches its maximum (see 
inset in FigjTJi) . r,-j > means that j turns before i, 
and vice versa. In absence of noise time ordering rela- 
tion -TOR - requires that = + Tkj , for each triplet 
i, j, k. We check robustness of this relation with respect 
to noise in all our flocks and find a relatively small spread 
of the data along the identity line (Fig. [4]). We recall that 
Tjj uses the full trajectory information of the two birds, 
which in turn is the product of the raw field data, of the 
tracking method and of the time-discrete data filtering. 
Hence, by proving that the TOR violation is low, we 
give a rather strong proof of reliability of our entire ex- 
perimental method. The quality of our TOR consistency 
test can be fully appreciated when we compare turning 
with non-turning flocks. If there is no turn is just a 
random number, so temporal consistency is strongly vi- 
olated and the TOR consistency plot really looks quite 
different (Fig. [4j lowest-right panel). 

Ranking. In absence of noise TOR would be exactly 
satisfied and ranking would be trivial. However, noise 
introduces some violations of TOR and we are in a simi- 
lar case as sport ranking, where player i may win over k, 
k may win over j, but i may lose to j, introducing some 
frustration. For every bird i, we say that i 'wins' over 
j if Tij < 0, in which case we set Jij = 1; conversely, i 
'loses' to j if Tij > 0, in which case we set J 



then define the total score of 



-1. We 

i as, (pi = J2j^i J ij- Given 
that we are in a round-robin tournament, it makes sense 
to rank the birds simply according to the scores (pi |49j . 
Thanks to the low violation of TOR, this score ranking 
already gives very small frustration (defined as the num- 
ber of cases in which i ranks higher than j, but i has lost 
to j). More refined rankings can be obtained by using 
probabilistic methods |49) . In our case, these methods 
decrease only marginally the (already low) frustration. 

Absolute turning time. The absolute turning time 
ti for each bird i is the delay with respect to the top bird 
in the ranking, i.e. the first to turn (f top = 0, r top = 1). 
However, to reduce the statistical error on tj introduced 
by TOR violations we define U using the mutual delay 
with respect to all birds j better ranked than i, 



+ Tij) 



n > l 



(M2) 



(a i (i)-a 3 (t-T)>-(a l (t))-(a J (i-T)) 



where (•) indicates a time average/integral and o~i is the 
fluctuation of a^(t) during the turn, 



va^m - (a a (t)> 2 ) . 



(Ml) 



Clearly, if there were no TOR violations, we would sim- 
ply have tj = T j,top- In the presence of noise, though, 
definition ( M2 ) is a more reliable estimate of U . By 
plotting Tj vs. ti for all birds in the flock, we obtain 
the ranking curve, r(t), which is reported for several of 
our flocks in Fig. [5] together with the propagation curve, 
x(t) = [r(t)/p] 1 '*. 
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Turning delay f (s) Time f (s) 

FIG. 5. Ranking and propagation. The ranking curve, r(t) and the propagation curve, x(t) = [r(i)/p] 1//3 , are reported for 
several turning flocks in our pool of data. 
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